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Abstract 

We propose to employ the quasiisotropic metamaterial (QIMM) slab to construct a polarization 
insensitive lens, in which both E- and H-polarized waves exhibit the same refocusing effect. For 
shallow incident angles, the QIMM slab will provide some degree of refocusing in the same manner 
as an isotropic negative index material. The refocusing effect allows us to introduce the ideas 
of paraxial beam focusing and phase compensation by the QIMM slab. On the basis of angular 
spectrum representation, a formalism describing paraxial beams propagating through a QIMM slab 
is presented. Because of the negative phase velocity in the QIMM slab, the inverse Gouy phase 
shift and the negative Rayleigh length of paraxial Gaussian beam are proposed. We find that the 
phase difference caused by the Gouy phase shift in vacuum can be compensated by that caused by 
the inverse Gouy phase shift in the QIMM slab. If certain matching conditions are satisfied, the 
intensity and phase distributions at object plane can be completely reconstructed at image plane. 
Our simulation results show that the superlensing effect with subwavelength image resolution could 
be achieved in the form of a QIMM slab. 
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I. INTRODUCTION 



About forty years ago, Veselago firstly introduced the concept of left-handed material 
(LHM) in which both the permittivity e and the permeability /x are negative l|. He predicted 
that LHM would have unique and potentially interesting properties, such as the negative 
refraction index, the reversed Doppler shift and the backward Cerenkov radiation. Veselago 
pointed out that electromagnetic waves incident on a planar interface between a regular 
material and a LHM will undergo negative refraction. Hence a LHM planar slab can act as 
a lens and focus waves from a point source. LHM did not receive much attention as it only 
existed in a conceptual form. After the first experimental observation of negative refraction 
using a metamaterial composed of split ring resonators (SRRs) , the study of such 
materials has received increasing attention over the last few years. While negative refraction 
is most easily visualized in an isotropic metamaterial negative refraction can 

also be realized in photonic crystals 0, H, B, and anisotropic metamaterials Q 
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19] have also been reported. 



Recently, Pendry extended Veslago's analysis and further predicted that a LHM slab 
can amplify evanescent waves and thus behaves like a perfect lens 20]. He proposed that 
the amplitudes of evanescent waves from a near-field object could be restored at its image. 
Therefore, the spatial resolution of the superlens can overcome the diffraction limit of con- 
ventional imaging systems and reach the subwavelength scale. The great research interests 
were initiated by the revolutionary concept. More recently, the anisotropic metamaterials 
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23j,|24j. Although 



have been proved to be good candidates for slab lens application 
the focusing is imperfect, the substantial field intensity enhancement can readily be ob- 
served. In these cases, the anisotropic metamaterials under consideration are characterized 
by a hyperboloid dispersion relation, and the focusing is restricted to either E- or H-polarized 
radiation. The recent development in quasiisotropic metamaterial (QIMM) offers us further 
opportunities to extend the previous work and further predict that both E- and H-polarized 
waves can be refocused. 

The main purpose of the present work is to construct a polarization insensitive lens by 
a QIMM slab. For shallow incident angles the QIMM slab will provide some degree of 
refocusing in the same manner as an isotropic LHM slab. We are particularly interested in 
exploiting polarization insensitive effect. First, starting from the representation of plane- 
wave angular spectrum, we derive the propagation of paraxial beams in the QIMM slab. 
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Our formalism permits us to introduce ideas for beam focusing and phase compensation 
of paraxial beams by using the QIMM slab. Next we want to introduce the inverse Gouy 
phase shift and negative Rayleigh length when waves propagating in the QIMM slab. As an 
example, we obtain the analytical description for a Gaussian beam propagating through a 
QIMM slab. We find that the phase difference caused by the Gouy phase shift in vacuum 
can be compensated by that caused by the inverse Gouy phase shift in the QIMM slab. 
If certain matching conditions are satisfied, the intensity and phase distributions at object 
plane can be completely reconstructed at the image plane. Finally, we will discuss what 
happen when the evanescent wave transmission through the QIMM slab. 



II. POLARIZATION INSENSITIVE METAMATERIAL 



Before we consider the polarization insensitive lens, we first analyze what is the QIMM. 
For anisotropic materials, one or both of the permittivity and permeability are second-rank 
tensors. In the following we assume that both the permittivity and permeability tensors are 
simultaneously diagonalizable: 
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where Ej and fij are the relative permittivity and permeability constants in the principal 
coordinate system (j = x,y,z). It should be noted that the real anisotropic metamaterial 



constructed by SRRs is highly dispersive, both in spatial sense and frequency sense 
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221 ] . So these relative values are functions of the angle frequency d. 



Following the standard procedure, we consider a monochromatic electromagnetic field 
E(r, t) = ife[E(r) exp(— iut)] and B(r,t) = ite[B(r) exp(— iut)] of angular frequency uj 
incident from vacuum into the anisotropic metamaterial. The field can be described by 



Maxwell's equations 



25) 



dB 

VxE = -aP 



B = fi fi • H, 

D = e £ ■ E. (2) 
The previous Maxwell's equations can be combined in a straightforward way to obtain the 
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FIG. 1: (Color online) (a) The QIMM with ellipsoid wave- vector surface; (b) The QIMM with 
double-sheeted wave-vector surface. 

well-known equation for the complex amplitude of the electric field, which reads 

Vx( ^' VxE) + ^ = ' (3) 

where c is the speed of light in vacuum. 

In the principal coordinate system, Maxwell's equations yield a scalar wave equation. In 
free space, the accompanying dispersion relation has the familiar form 

k x + k y + k z = — , (4) 

where kj is the j component of the incident wave vector. 

We note that the Maxwell's equations are symmetrical in electric and magnetic fields. So 
as far as Maxwell's equations is concerned, what we can do for electricity we can also do for 
magnetism. To achieve the polarization insensitive effect, we will focus our interesting on 
the anisotropic metamaterial, in which the permittivity and permeability tensor elements 
satisfy the condition: 

E^l EM. EM C (C> 0), (5) 

fJL x (u) HyiU)) H Z (UJ) 

where C is a constant. A careful calculation of the Maxwell's equations gives the dispersion 
relation: 



_£L + jL + jL_^Vjg_ + jL + jL_^' 1 = o, (6 ) 

E-yl^z ^■xP'z £y[^x C J \E z fAy Ezf^x ^-xf^y C 
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where qj represents the j component of transmitted wave-vector. The above equation can 
be represented by a three-dimensional surface in wave-vector space. This surface is known 
as the normal surface and consists of two shells 25j . Under the condition of Eq. (J5J) , we can 



find E- and H-polarized waves exhibit the same wave-vector surface. Thus the anisotropic 
medium also be regard as QIMM [26 



[27] . Clearly, we can find the dispersion surface has 
the following two types: ellipsoid and double-sheeted hyperboloid, as show in Fig. [TJ 

We are currently investigating the possibilities for the manufacture of the QIMM. In 
fact, it is now conceivable that a metamaterial can be constructed whose permittivity and 
permeability values may be designed to vary independently and arbitrarily throughout a 
metamaterial, taking positive or negative values as desired. Hence the permittivity and 
permeability tensor elements of QIMM can be controlled by modulating the length scale 
of the SRRs. The permittivity £j(uj) and permeability /J,j(u) can be approximated by the 
Lorentz model. We trust that the QIMM can be constructed, because similar technology 
has been exploited in anisotropic metamaterial , ll7J, [19] . In addition, photonic crystals 



might be a good candidate for constructing a QIMM. The periodicity in photonic crystals 
is on the order of the wavelength, so that the distinction between refraction and diffraction 
is blurred. Nevertheless, many novel dispersion relationships can be realized in photonic 
crystals, including ranges where the frequency disperses negatively with wave vector as 
required for a negative refraction {7, [§], O, lo| . 

Now we want to enquire: whether E- and H-polarized exhibit the same propagation 
characteristic. To answer the question we first discuss the transmission of wave vector. We 
choose the z axis to be normal to the interface, the x and y axes locate at in the plane of the 
interface. The ^-component of the transmitted wave vector can be found by the solution of 
Eq. (jSJ), which yields 



q E z = a x je y fi x k 2 - e y fi x ( + ) , (7) 



Qz = v\le x Hyk% - e x fi y ( + , (8) 

for E- and H-polarized waves, respectively. Here ko = uo/c is the wave number in vacuum 
and a = ±1. This choice of sign ensures that power propagates away from the boundary to 
the +z direction. 

Without loss of generality, we assume the wave vector locate at the x — z plane (k y = 
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q y = 0). The incident angle of light is given by 



9j = tan 



(9) 



The values of refractive wave vector can be found by using boundary conditions and disper- 
sion relations. The refractive angle of the transmitted wave vector or phase of E- and H- 
polarized waves can be written as 



P. 



tan 



of 



f3p = tan 



9: 



11 



(10) 



Substituting Eqs. (j7j) and (jSj) into Eq. ffTUl) . we can easily find that E- and H- polarized 
waves propagates with the same wave vector or phase velocity. 

In the next step, let us discuss the transmission of energy flux. It should be noted the 
actual direction of light is determined by the time-averaged Poynting vector S = |Re(E* x 
H). For E- and H-polarized waves, the transmitted Poynting vector St is given by 



Re 



Re 



T E E 0lx , T E E l^ 
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1 H E oQ: 



TZ T?Z r ,l 



(11) 



(12) 



2ue z 2a;e 3 

where T% and Tjj are the transmission coefficients for E- and H-polarized waves, respectively. 
The refraction angle of Poynting vector of E- and H- polarized incident waves can be obtained 

as 



0, 



tan 



qE 



Tx 



qE 



/3§ = tan" 



ch n 

CH 



(13) 



Combining Eqs. (flOj) and (Tl3|) we can easily find that E- and H- polarized waves have the 
same Poynting vector. As for the QIMM slab, the refraction at the second interface can be 
investigated by the similar procedures. 

By now, we know that E- and H- polarized waves propagate with same wave vector 
and Poynting vector. It is a significantly different property from general anisotropic media. 
Note that there is a bending angle between q and S, and therefore q, E and H do not 
form a strictly right-handed or left-handed system in QIMM. Hence it is also different 
from isotropic media. For this reason, this kind of special anisotropic media is regarded as 
quasiisotropic. It should be mentioned that if C < in Eq. (JSJ), E- and H-polarized waves 
will exhibit the same single-sheeted dispersion relation. While the two polarized waves will 
undergo different amphoteric refraction, the special anisotropic media cannot be regarded 



as quasiisotropic 



27] . 




FIG. 2: The frequency contours of isotropic and quasiisotropic media: (a) The circle and ellipse 
represent the frequency contours of vacuum and quasiisotropic media, respectively. Both the wave 
vector and the energy flow exhibit negative refraction, (b) The circle and hyperbola denote the 
frequency contours of vacuum and quasiisotropic media, respectively. The wave vector undergoes 
a positive refraction, while the energy flow undergoes a negative refraction. 

Now we are in the position to study the negative refraction in the QIMM. Unlike in 
isotropic media, the Poynting vector in The QIMM is neither parallel nor antiparallel to 
the wave vector, but rather makes either an acute or an obtuse angle with respect to the 
wave vector. In general, to distinguish the positive and negative refraction in QIMM, we 
must calculate the direction of the Poynting vector with respect to the wave vector. Positive 
refraction means q x • St > 0, and negative refraction means • St < 0. From Eqs. ( TTTT) 
and (JI21) we get 

rn2 7712 ^2 rp2 Tj2„2 

The negative refraction phenomenon is one of the most interesting properties of the QIMM. 
We can see that the refracted waves will be determined by [i z for E-polarized incident waves 
and e z for H-polarized incident waves. 

Because of the importance of negative refraction in refocusing effect, we are interested 
in the two types of QIMM, which can formed from appropriate combinations of material 
parameter tensor elements. 

Type I. In this case all of the Sj and fij are negative. The frequency contour is an ellipse as 
shown in Fig. [2(a). Here k 2 ■ q z < and q x • St < 0, so the refraction angle of wave vector 
and Poynting vector are always negative. 
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FIG. 3: The ray tracing picture showing the focussing by QIMM slab. The QIMM slab is sur- 
rounded by vacuum in region 1 and region 3. The solid line and dash-dotted lines are the theoretical 
objective and focusing planes, respectively. 

Type II. In this case e x > 0, e y > and e z < 0. The frequency contour is a double-sheeted 
hyperbola as depicted in Fig. 12(b). Here k z • q z > and q x ■ St < 0. It yields that the 
refraction of Poynting vector refraction is always negative even if the wave- vector refraction 
is positive. 

As noted above, the Poynting vector will exhibit negative refraction in the two types of 
QIMM. The negative refraction is the important effect responsible for the slab lens. Hence, 
the two kind of QIMM can be employed to construct a polarized insensitive lens. 



III. THE PARAXIAL MODEL OF BEAM PROPAGATION 

In this section, we consider the slab lens constructed by the QIMM. As depicted in Fig. |3j 
the QIMM slab in region 2 is surrounded by vacuum in region 1 and region 3. A point source 
is place the object plane z — 0. A single ray will pass the interfaces z = a and z = a + d 
before it reaches the image plane z = a + b + d. Let us investigate what happens when 
the single ray pass through the QIMM slab. Because of the anisotropy the ray incident for 
different angle will exhibit different image plane. 

First, we explore the aberration effect caused by the anisotropic effect. Subsequent cal- 



S 



dilations of Eq. (fl3j) give the relationship between sin 9j and sin (5 S by 

sin 2 /3f= - A^L ; (15 ) 



9 „ ff £~ sin 2 #r , , 

Sil1 = 5 f \ ■ 2n ■ 16 

Note that the expressions are sightly different from those in conventional uniaxial crystal 24]. 
It is instructive to compare these results with Snell's law, which for refraction from vacuum 
into an isotropic medium, gives sin (9// sin /3s = n, where n represents the refractive index of 
the refracting medium. Evidently, in the present case we can find the relationship between 
sin 9 1 and sin/?,, is nonlinear, which is caused by the anisotropic effect e x ^ e z or \h x ^ \i z . 
The nonlinear relationship will result in a significant aberration effect in the image plane, 
hence the achievable resolution of image is limited. 

Next, we want to discuss another aberration effect caused by frequency dispersion. For 
a certain incident angle, the ray with different frequency will exhibit different image plane 
as shown in Fig. HI We assume that both e z (u) and [i z (u) be approximated by the Lorentz 
model. The other tensor elements components are approximated as constants. The resonate 
frequency, plasma frequency, and damping constant are identical in all respects to those 
utilized in Ref. [lTj] - Ignoring the metallic structure, the other tensor elements assume 
the values of the background material which is dominantly air. Obviously the frequency 
dispersion of e z {uS) and fi z (u) will place some practical limitations on the resolution of 
image. Fortunately the limitations can be reduced, if the ray incident at a small angle. 
The significant effect has be illustrated in Fig. HI It is interesting to noted that the image 
distance will slowly vary with frequency in a certain band. 

Remarkable as the QIMM slab lens is it suffers from some problem: how to cancel the 
aberration effect? Does the polarization insensitive lens have any use, if the image is not 
imperfect? The main effect of anisotropy associated with QIMM will limit the resolution 
of the image. Although the image is imperfect, we trust that the reconstruct effect of 
intensity and phase in paraxial regime will lead to some potential applications. Such as 
the QIMM slab can be used to provide phase compensation and beam focusing in cavity 
resonator. Furthermore, sightly anisotropy in QIMM slab lens can improve the beam param- 
eter. Most importantly, the QIMM can be used to design polarization-insensitive modulators 
and polarization-insensitive all-optical switching in fiber communication system. Hence it is 
very desirable to investigate the polarization insensitive effect in paraxial regime. 
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FIG. 4: For a certain incident angle, the ray with different frequency / = lo/2tt will exhibit different 
image plane. Note that for smaller incident angles, the QIMM slab will provide some degree of 
frequency insensitive in a certain GHz band. 

le field can 
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From a mathematical point of view, the approximate paraxial expression for t 
be obtained by the expansion of the square root of q z to the first order in |q±|/? 
which yields 

q E = a^/TJTko - ( a ^ £y ^ x ^ x _|_ a \ /£ y^ x K \ ^yj) 
V 2£yV' Z ko 2e x [i z k J 

Qz - a^ko - {-^^ + ^J-) . (18) 



where we have introduced the boundary condition = k^. From Eqs. (IT7I) and ([18]) we 
can easily find that for shallow incident angles the QIMM slab will provide some degree 
of refocusing in the same manner as an isotropic LHM. Hence the aberration effect can be 
cancelled in paraxial beam region. The interesting property allow us to introduce the idea 
to construct a QIMM slab lens in paraxial beam region. 

Equation can be conveniently solved by employing the Fourier transformations, so the 
complex amplitudes in QIMM for E- and H-polarized beams can be conveniently expressed 
as 

E(r ± ,z)= J d 2 k ± E(k ± )exp\ik ± -v ± + iqfz]. (19) 
H(r_i_, z) = J d 2 k ± H(k ± ) exp[ik ± ■ r ± + iqfz). (20) 
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Here = xe x + ye y , kj_ = k x e x + k y e y , and ej is the unit vector in the j-direction. 
Substituting Eqs. (fT7|) and f|T8l) into Eq. (fT9l) and fl20|) . respectively, we obtain 



E(r ± , z) = ex.p(icry/e x [iyk z) J d 2 k A 



x exp 



-t ± -( °i^kl + ^^k 2 
\2e y /j, z k 2e x fj, z k y 



E(kj 



(21) 



H(r ± ,z) 



ex.p(iay/e x fi y k z) J d 2 k 
ik ± • r ± - 



x exp 



-hi + 



Hfk 



(22) 



2e z /i y ko" x ' 2e z [i x k v y 

The fields E(k±) and H (kj_) in Eqs. (1211) and (1221) are related to the boundary distributions 
of the fields by means of the relation 



Efk, 



Hfk, 



Gpr_LE(r_L, 0) exp[zk ± • r J 



Gpr_LH(r_L, 0) exp[ik_L • r J 



(23) 



(24) 



for E- and H-polarized beams, respectively. Evidently, Eqs. (1231) and (1241) are standard two- 



dimensional Fourier transform [28| . In fact, after the field distribution in the plane z = 
is known, Eqs. (T2T1) and (l2"2"j) provide the expression of the E- and H-polarized field in the 
space z > 0, respectively. 

Since our attention will be focused on beam propagating along the +z direction, we can 
write the paraxial fields as 



E(r ± , z) = A e (t±, z) exp(iay/e y fx x koz) 



(25) 



H(r_L, z) = A H (rj_, z) exp(ia y /e x ~jl y ~k z), (26) 
where the field A(r±, z) is the slowly varying envelope amplitude which satisfies the parabolic 



equation: 



\ d ( cr^/eyji^ d 2 + <J<Je y ~jI x ~ d 2 
dz \2e y /j, z k dx 2 2e x /i z k dy 7 



'. d ( cr^/e^ d 2 cr^/e^ d 2 
dz \2e z /j, y k dx 2 2e z /j, x k dy 2 



A E (r x ,z) = 0, 



A H (r ± ,z) = 0. 



(27) 



(28) 
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Under the quasiisotropic condition of Eq. (jSj), we can easily find that E- and H-polarized 
paraxial filed exhibit the same propagating characteristics in paraxial regime. The interest- 
ing properties allow us to introduce the idea to construct a polarization lens by QIMM slab. 
For simplify, we introduce the effective refraction indexes: 

n x = a n y = a (29) 

From Eqs. fl271) and fl28l) we can find that the field of paraxial beams in QIMM can be written 
in the similar way to that in regular material, while the sign of the effective refraction index 
could be reverse. To simplify the proceeding analyses, we will focus our attention on the 
QIMM with ellipsoid frequency contour. 



IV. BEAM FOCUSING BY POLARIZATION INSENSITIVE LENS 



In the previous section we have understood both E- and H-polarized beams have the 
same propagation characteristic in QIMM slab. Hence we do not wish to get involved in 
the trouble to discuss the focusing effect of two polarized waves. Instead, we will investigate 
the analytical description for E-polarized beam with a boundary Gaussian distribution. 
This example allows us to describe the refocusing features of beam propagation in QIMM 
slab. To be uniform throughout the following analysis, we introduce different coordinate 
transformations z*(i = 1,2,3) in the three regions, respectively. First we want to explore 
the field in region 1. Without any loss of generality, we assume that the input waist locates 
at the object plane z = 0. The fundamental Gaussian spectrum distribution can be written 
in the form 



Ei(kj 



WqEq 

V2^ 



exp 



(30) 



where w is the spot size. The Rayleigh lengths give by Zr = k WQ/2. By substituting 
Eq. (!30|) into Eq. (fl9l) . the field in the region 1 can be written as 



Ei(r ± ,^) 



WqEq 



exp 



X 



+ 



w 



hi 



+ #i 



(31) 



wi x =w \ 1 + 



Wly = W \ 1 + 



J }y_ 



(32) 



Here we have chosen different waists, W\ x and Wi y , in order to deal with a more general 



situation. Because of the isotropy in vacuum, we can easily obtain z\ 



z and 



Wi x = Wi v . The corresponding Rayleigh lengths give by L 



lx 



ny — Z R- 
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We are now in a position to investigate the field in region 2. In fact, the field in the first 
boundary can be easily obtained from Eq. fl3Tl) by choosing z = a. Substituting the field 
into Eq. ([23]) . the angular spectrum distribution can be obtained as 



E 2 (kj 



WqEq 

V2^ 



exp 



(33) 



For simplicity, we assume that the wave propagate through the boundary without reflection. 
Substituting Eq. ( |33l) into Eq. ( f2~lT) . the field in the QIMM slab can be written as 



E 2 (r ± ,4 



WqEq 



v /W 2x W 2 y 



exp 



x 



+ 



w 



2y 



+ iip2 



(34) 



W 2a 



W \ 1 + 



Z 2x 



J 2x 



, W2y = Wo\ll + 



3l 



(35) 



Here z 



2x 



'1 — n x )a and = z 



'\ — n y )a. The interesting point we want to stress 
is that there are two different Rayleigh lengths, L 2x = n x k wl/2 and L 2y = n y k WQ/2, 
that characterize the spreading of the beam in the direction of x and y axes, respectively. 
A further important point should be noted that we have introduce the negative Rayleigh 
length. The inherent physics underlying the negative Rayleigh length is the waves undergo 
a negative phase velocity in the QIMM slab. As can be seen in the following, the negative 
Rayleigh length will give rise to the corresponding reverse Gouy phase shift. 

Finally we want to explore the field in region 3. The field in the second boundary can be 
easily obtained from Eq. (l34|) under choosing z = a + d. Substituting the field into Eq. ( !23f) . 
the angular spectrum distribution can be written as 



E,(k, 



WqEq 



exp 



n x koWQ + 2in x a + 2id 2 n y k WQ + 2in y a + 2id 2 

; z k„ ~\- ■ ■ k„, 



An x k 



An y k 



Substituting Eq. ( 1361) into Eq. ( 12T1) . the field in the region 3 is given by 



(36) 



E3(ri,*5) 



WqEq 



^w 3x w 3y 



exp 



x 



+ 



+ #3 



(37) 



w 3x = w \ 1 + 



J 3x 



W 3 y = Wq\ 1 + 



J3y 



(38) 



Here z\ x = z — (1 — l/n x )d and z 3y = z — (1 — \/n y )d. The corresponding Rayleigh lengths 



give by L 



3x 



J 3y 



kQW^/2, that denote the beam exhibit the same diffraction distance in 
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(a) (b) 




FIG. 5: (Color online) The numerically computed intensity distribution in object and image planes, 
(a) the intensity distribution for normal Gaussian beam in object plane. The in tensity distribution 
in image plane for Gaussian beam propagating through the QIMM slab with different anisotropic 
parameters: (b) n x = —1, n y = —2. (c) n x = —2, n y = —1. (d) n x = —1, n y = —1. We can easily 
find the intensity distribution at the object plane can be completely reconstructed at the image 
plane. 

the direction of x and y axes. The effect of the anisotropic diffraction is that these two beam 
widths keep their difference even if the Rayleigh lengths, L 3x and L 3y , are equal, implying 
that generally the Gaussian beam is astigmatic. 

Up to now, the fields are determined explicitly in the three regions. Comparison of 
Eq. (1341) . Eq. (1371) with Eq. (I3T1) shows that the field distributions in region 2 and region 3 
may no longer remain Gaussian. We take the image position z = a + d + b to be the place of 
the second focusing waist. For the purpose of illustration, the intensity distribution in object 
plane is plotted in Fig. [5]^a). In general, the shape of intensity distribution is distorted in 
image plane as shown in Fig. [5](b) and Fig. [5](c). Careful evaluation of Eq. fl371) reveals that 
the secret underlying the intensity distortion is the anisotropic diffraction. 

Now, the most obvious question is whether the intensity distribution at the object plane 
can be completely reconstructed at the image plane. In the next step, we want to explore 
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the matching condition of focusing. We can easily obtain the place of the focusing waist by 
choosing z* = 0. Let us assume the incident beam waist locates at plane z = 0. To eliminate 
the astigmatic effect, the beam waists should locate at the same place, namely z\ x = z£ y . 
Using these criterions, the matching condition for focusing can be written as 

e y \i z {a + b) + a<Je x /j y d = 0, e x = e y . (39) 

Under the focusing matching condition, the intensity distribution at the object plane can be 
completely reconstructed at the image plane as shown in Fig. [5](d) . A further point should be 
noted is that the thickness of the QIMM slab should satisfy the relation d > ae y fi z a/ \ye x fi y , 
otherwise there is neither an internal nor an external focus. 



V. PHASE COMPENSATION BY POLARIZATION INSENSITIVE LENS 



In this section, we attempt to investigate the matching condition for phase compensation. 
In isotropic LHM, plane waves can propagate with negative phase velocity directed opposite 
to the direction of Poynting vector. Hence the phase difference can be compensated by 



the LHM slab 



271 ]. However, the negative tensor parameters associated with QIMM 
provides a wealth of opportunities for observing and exploiting negative phase-velocity be- 
havior. 

First let us investigate the phase distribution in region 1. A more rigorous calculation of 
Eq. (131]) gives 

* = k °> + (%r + Hr)-*- (40) 

Rix = 4 X + % Riv = 4 y + (41) 

Z lx Z ly 



1 / Z* Z* \ 

$i = — I arctan — — + arctan — — I . 

2 V Li x Li y J 



(42) 



Here R\ x and R\ y are the radius of curvature. Because of the isotropy in vacuum, we can 
easily find R(z* x ) = R(z* ). The Gouy phase shift in vacuum is given by $i. 

Next, we attempt to explore the phase distribution in region 2. Matching the boundary 
condition, the phase term in Eq. (13"4"|) can be written as 

ip 2 = ha + a^/EynMz -a) + — — + — $ 2 , (43) 
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FIG. 6: The phase difference caused by the Gouy phase shift in vacuum can be compensated by 
that caused by the inverse Gouy phase shift in the QIMM slab. The phase fronts of Gaussian beam 
(solid lines) differ from those of a perfect spherical wave (dashed lines). 

R2 X = 4x + |K R2y = 4y + (44) 

Z 2x Z 2y 



2 



1 / Z* Z* \ 

-- I arctan — + arctan . (45) 

2 V L,2 X L<2y J 



The Gouy phase shift in QIMM is given by Eq. (j4"5l) . We should mention that there are 
two different radius of curvature, R2 X and R2 y , that characterize the beam undergo different 
diffraction effects in the direction of x and y axes, respectively. 

Now, we are in the position to explore the phase distribution in region 3. Analogously, 
we make some serious calculation of Eq. (|37|) . then obtain the phase distribution 

^3 = (Ty/e y fi x k d + k (z -d) + + - $ 3 , (46) 
Rs* = 4x + |k ifey = 4y + ^ (47) 

z 3x z 3y 



1 / z* z* \ 

— I arctan + arctan . (48) 

2 V Ls x Lz y J 



The radius of curvatures are given by Eq. f|471) and the corresponding Gouy phase shift is 
given by Eq. (1481) . The anisotropic effect result in the two radius of curvatures keep their 
difference even if the Rayleigh lengths are equal. 
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(a) m 




FIG. 7: (Color online) The numerically computed phase distribution in object and image planes, 
(a) The phase distribution in object plane. The phase distribution in image plane after the Gaussian 
beam propagating through the QIMM slab with different anisotropic parameters: (b) n x = — 1, 
n y = —2. (c) n x = —2, n y = —1. (d) n x = —1, n y = —1. The phase distribution can be completely 
reconstructed at the image plane. The parameters are the same as in Fig. [5j 

It is known that an electromagnetic beam propagating through a focus experiences an 
additional n phase shift with respect to a plane wave. This phase anomaly was discovered 



by Gouy in 1890 and has since been referred to as the Gouy phase shift [32|, |33J]. It should be 



mentioned that there exists an effect of accumulated Gouy phase shift when a beam passing 



through an optical system with positive index [34], |35|, [36j| . While in the QIMM slab system 
we expect that the phase difference caused by Gouy phase shift can be compensated by that 
caused by the inverse Gouy shift in the QIMM slab. 

We might suspect whether the phase difference caused by the Gouy phase shift in vacuum 
can be compensated by that caused by the inverse Gouy phase shift in QIMM slab. To obtain 
the better physical picture, the schematic distribution of phase fronts are plotted in Fig. El 
The phase fronts of a focused Gaussian beam are plotted with solid lines, and the phase 
fronts of a perfect spherical wave are depicted with the dashed lines. The phase difference 
on the optical axis is caused by the Gouy phase shift. The inherent secret underlying the 
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reverse Gouy phase shift in the QIMM slab is the waves undergo a negative phase velocity. 

Let us investigate what happens if we consider the phase difference caused by the Gouy 

shift. Under the focusing matching conditions, the phase difference caused by the Gouy 

phase shift in the three regions are 

. _ a 
Aq>! = — arctan — , 

. a b 
A<P 2 = arctan h arctan — , 

Zr Zr 

A$ 3 = -arctan — . (49) 

zr 

The first and third equations dictate the phase difference caused by the Gouy shift in regions 
1 and 3, respectively. The second equation denotes the phase difference caused by the inverse 
Gouy phase shift in the QIMM slab. Subsequent calculations of Eq. f l4"9"j) show 

A$i + A$ 2 + A$ 3 = 0. (50) 

This implies that the phase difference caused by the Gouy phase shift in vacuum can be 
compensated by the counterpart caused by the inverse Gouy phase shift in QIMM slab. 
Therefore the condition for phase compensation can be simply written as 

(a + b)k + <jy/e y {i x dk = 0. (51) 

The first term in Eq. ( 15T1) is the phase deference caused by the plane wave in vacuum, and 
the other term is the phase deference caused by the plane wave in the QIMM slab. 

For the purpose of illustration, the phase distribution in object plane is plotted in 
Fig.[7Ja). Generally, the phase distributions in image plane is distorted as shown in Fig. [7Jb) 
and Fig. [7(c). As mentioned above, the phase distortion is caused by the effect of anisotropic 
diffraction. To cancel the phase distortion, the beam waists should locate at the same place, 
namely z\ x = z\ y . Under the phase matching condition, the phase distribution at the object 
plane can also be completely reconstructed at the image plane as depicted in Fig. [7J(d). 

Now an interesting question naturally arises: whether the matching conditions of focusing 
and the phase compensation can be satisfied simultaneously. Clearly, if we seek a solution 
satisfying Eqs. fl39l) and fl5TT) . the only possibility is 

&x = £yi ^yf^z = 1- (52) 

Under the matching conditions, the intensity and phase distributions at the object plane 
can be completely reconstructed at the image plane. 
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It should be mentioned that, for the QIMM slab with double-sheeted hyperboloid wave- 
vector surface, both E- and H-polarized beams can also exhibit the same intensity and 
phase reconstructed effect. Because of the positive phase velocity embedded in this type of 
QIMM, the paraxial beam will experience the positive Rayleigh distance and Gouy phase 
shift. Therefore the accumulated phase delay effect give rise to a large phase deference 
between the object and image planes. 



VI. THE TRANSMISSION OF EVANESCENT WAVES 

In this section, we discuss under what conditions anomalous transmission will occur 
when an evanescent wave transmitted through the QIMM slab. It is well known that when 
an evanescent wave is transmitted through a slab of regular media with simultaneously 
positive permittivity and permeability, the amplitude of the transmitted wave will decay 
exponentially as the thickness of the slab increases. While an evanescent wave transmitted 
through an isotropic LHM slab, the amplitude of the transmitted wave would be amplified 
exponentially. This anomalous transmission of evanescent waves is a very peculiar property 
of LHM and it may lead to subwavelength image {20I, [i?! ■ 

Now we will explore what happen when the evanescent wave transmission through the 
QIMM slab. For E-polarized incident waves, the incident and reflected fields in region 1 can 
be written as 

Ei = E e y exp[i(k x x + k z z)\ + REE e y exp[i(k x x — k z z)}, (53) 

where Re is the reflection coefficient. Some of the incident wave is transmitted into the 
QIMM slab, and conversely a wave inside the QIMM slab incident on its interfaces with the 
surrounding vacuum also experiences transmission and reflection, so the electric field of the 
wave inside the slab is given by 

E 2 = rE e y exp[i(q^x + qfz)] + tE e y exp[i(qfx - qfz)]. (54) 



Here r and t are coefficients which need to be determined by boundary conditions [12 ]. 
Matching the boundary conditions for each wave- vector component at the plane z = a + d 
gives the propagation field in the form 

E 3 = T E E e y exp[iqfx + qf{z - d)}, (55) 
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FIG. 8: An evanescent wave (k x = 2ko) interacting with a QIMM slab (a) d = 0.2A (b) d = 0.4A. 
The QIMM with different values of absorption 7 e = j m = 7, ranging from 0.001 to 0.1. The QIMM 
slab is supposed to be well satisfied the conditions: e x = e y = —0.8 — ji and fi z = —1.25 — ji. The 
vertical gray lines denote the two surfaces of the QIMM slab. 

where Tg is the overall transmission coefficient. The ^-component of the wave vectors of 
evanescent waves can be found by the solution of Eq. ([6]), which yields 



qf = iJ—q 2 x - e y fi x kl qf = iJ j-q* - e x ii y k%, (56) 
V /Az V z 

for E- and H-polarized waves, respectively. 

By matching the electric and magnetic fields at the two interfaces between the QIMM slab 
and the surrounding vacuum, the coefficients in Eqs. ( |53l) . ( 1541) . and (|55l) can be determined. 
We can get that the overall transmission through both surfaces of the QIMM slab is given 
by 



E I ■■ ^ ' ~,E\2 _ (r,E _ ,, U \2 On r,EA\-> 



4/J> x k z qf exp(iqfd) 
i^xkz + qf ) 2 ~ (gf - fixkz) 2 exp{2iqfd) ' 
From Eq. fl57|) we can see that in general cases, when an evanescent wave is transmitted 

through a QIMM slab, its amplitude will decay exponentially as the thickness of the slab 

increases. But if the following conditions are satisfied: 

fi x < 0, fi z < 0, e y < 0, (58) 

e x ji z = 1, Eyfiz = 1, (59) 

the amplitude of the transmitted evanescent wave will be amplified exponentially by the 
transmission process through the QIMM slab. From Eq. (!57|) . we can see that if the con- 
ditions fl58l) and fl59l are satisfied, the overall transmission coefficient T E will be equal to 
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exp(\qf\d), hence the amplitude of the transmitted evanescent wave will be amplified expo- 
nentially as the thickness of the QIMM slab increases. Note that the conditions mentioned 



above are sightly different from those as Hu and Chui have obtained [121 ]. 

As for H-polarized evanescent waves, the overall transmission through both interfaces of 
the QIMM slab can be obtained by similar procedures, and we can get the overall transmis- 
sion coefficient: 

rp = 4e x k z q* exp(iq*d) 

From Eq. (1601) . we can see that in general cases, the amplitude of the transmitted H-polarized 
evanescent waves will also decay exponentially as the thickness of the QIMM slab increases. 
But if the following conditions are satisfied, the overall transmission coefficient Th will be 
equal to exp(\q^\d), and the amplitude of the transmitted H-polarized evanescent wave will 
be amplified exponentially by the transmission process through the QIMM slab. So for 
H-polarized evanescent waves, 

e x < 0, e z < 0, s y < 0, (61) 

e z fi x = 1, e z \i y = 1. (62) 
if conditions flBTl) and (l62l) are satisfied, the QIMM slab will enhance exponentially the 



transmitted waves. Comparing (1581) and (1591) with ( ISTT) and ( 1621) . we find that in the presence 
of QIMM, the conditions for E-polarized means automatically the conditions for H-polarized 
are satisfied. Hence the QIMM also exhibit the significant insensitive effect for evanescent 
waves. We stress that, for the QIMM slab with double-sheeted hyperboloid dispersion 
relation, the evanescent waves cannot be amplified. 

In the above analysis the permittivity and permeability tensor elements are assumed 
to be lossless. However the effect of absorption, necessarily present in such materials, may 



drastically suppress any evanescent amplifying wave into a decaying one [38|, |39( . Comparing 
Fig. [8]^a) with Fig. [8](b) suggests that large d is not favored for the amplification of the 
evanescent wave inside the QIMM slab. It is also clearly seen from Figs. [8](a) and^b) that 
a low reflection occurs at the interfaces, since the QIMM slab is supposed to be well satisfied 
the conditions ( l58|) and ( |59l) . So far we have shown the simulation results of the evanescent 
wave interacting with QIMM slabs with different absorptions and thicknesses. We find 
that the suppression of evanescent- wave amplification can be effectively relaxed by reducing 
the thickness of the QIMM slab. The numerical examples provide direct evidence that an 
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evanescent wave could be amplified in a QIMM slab with finite absorption. Unfortunately, 



the loss of a realistic metamaterial could not be reduced to a very small level [40|. Therefore, 
the thickness of the QIMM slab should be much smaller than the wavelength in order to 
realize the subwavelength imaging for such a system. 

As a result of the amplification of evanescent waves inside the lossy QIMM slab, super- 
lensing effect with subwavelength image resolution could be achieved practically. While the 
main effect of anisotropy associated with QIMM will limit the resolution of the image. Al- 
though the image is imperfect, we trust that the reconstruct effect of intensity and phase in 
paraxial regime will lead to some potential applications. Several recent developments make 
the polarization insensitive lens a practical possibility. Some time ago it was shown that a 



double-periodic array of pairs of para 
permeability in the optical range 41 



lei gold nanorods will exhibit negative permittivity and 



421 ] . Another extremely promising material has been 



previously explored in certain designs of photonic crystals, which can be effectively modelled 
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451 ] . Experimentally, the goal 



with anisotropic permittivity and permeability tensors 
to realize the insensitive lens, chiefly lies in reducing the loss of QIMM. Practical polarization 
insensitive lens will require the frequency independent, a great challenge to the designers is 
to realize the negative material parameters in a wide band. The recent developments lead 
us to be optimistic that the polarization insensitive lens can be designed in future. 



VII. CONCLUSIONS 



In conclusion, we have proposed how to employ the QIMM slab to create a polarization 
insensitive lens, in which both E- and H-polarized waves exhibit the same refocusing effect. 
For shallow incident angles the QIMM slab will provide some degree of refocusing in the 
same manner as an isotropic negative index material. We have investigated the focusing 
and phase compensation of paraxial beams by the QIMM slab. We have introduced the 
concepts of inverse Gouy phase shift and negative Rayleigh length of paraxial beams in 
QIMM. We have shown that the phase difference caused by the Gouy phase shift in vacuum 
can be compensated by that caused by the inverse Gouy phase shift in the QIMM slab. 
If certain matching conditions are satisfied, the intensity and phase distributions at object 
plane can be completely reconstructed at the image plane. The QIMM slab exhibits the 
significant insensitive effect for both transmitted and evanescent waves. Our simulation 
results show that the super lensing effect with subwavelength image resolution could be 



22 



achieved in the form of a QIMM slab. We wish the essential physics described in this paper 
will provide reference in the road to construct the polarization insensitive lens. We trust 
that the significant insensitive properties will lead to further novel effects and applications. 
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